Fock Space Representations for Non-Hermitian Hamiltonians 



David B. Fairlie 1 and Jean Nuyts 2 

David.Fairlie@durham.ac.uk, Jean.Nuyts@umh.ac.be 

1 Department of Mathematical Sciences, University of Durham 
South Road, DH1 3LE, Durham, England 

2 Physique Theorique et Mathematique, Universite de 

Mons-Hainaut 
20 Place du Pare, 7000 Mons, Belgium 



Abstract 

The requirement of Hermiticity of a Quantum Mechanical Hamiltonian, for the 
description of physical processes with real eigenvalues which has been challenged 
notably by Carl Bender, is examined for the case of a Fock space Hamilitonian 
which is bilinear in two creation and destruction operators. An interpretation of 
this model as a Schrodinger operator leads to an identification of the Hermitian 
form of the Hamiltonian as the Landau model of a charged particle in a plane, 
interacting with a constant magnetic field at right angles to the plane. When the 
parameters of the Hamiltonian are suitably adjusted to make it non-Hermitian, 
the model represents two harmonic oscillators at right angles interacting with a 
constant magnetic field in the third direction, but with a pure imaginary 
coupling, and real energy eigenvalues. It is now VT symmetric. Multiparticle 
states are investigated. 



1 Introduction 



The familiar Hamiltonians of Quantum Mechanics may be analysed for sym- 
metries either in terms of pure matrix algebra, or else in terms of a Fock space 
representation, which generally leads to a more physical interpretation of the 
mathematical manipulations. In particular, any Hamiltonian constructed 
from a Fock space of n fermionic creation and annihilation operators may 
be transcribed in terms of a finite dimensional matrix of dimension 2 2 x 2 a 
on account of the fact that such operators can always be constructed from 
gamma matrices, which admit a well known matrix representation. The in- 
verse of this construction, to obtain a Fock space representation of any finite 
dimensional Hermitian matrix can be done by embedding the matrix in a 
larger one of suitable dimension to admit representations of the canonical 
anti-commutation relations. In fact, in view of the existence of a matrix 
representation for Hamiltonians the transcription is merely an exercise in 
matrix algebra for undergraduates. However, the situation is more compli- 
cated when the Hamiltonian is no longer Hermitian. It is a central tenet of 
Quantum Theory and Quantum Field Theory that the Hamiltonian should 
be Hermitian, but this is not in fact necessary to maintain real eigenvalues, as 
has been demonstrated in recent years, notably by Bender and Boettcher.0. 
However, if the eigenvalues are all discrete and real and their eigenvectors 
span the full space, there will exist, in general, an infinity of Hermitian 
Hamiltonians (possibly infinite dimensional) which are unitarily equivalent 
to the diagonal (real) form. They are equivalent up to a more general (not 
necessarily unitary) transformation to the initial Hamiltonian. 

If all the eigenvalues are real and their eigenvectors span a subspace S of 
the initial space (say finite dimensional to make things simpler) there exists, 
in general, an infinity of Hamiltonians unitarily equivalent to a Jordan form 
(or a generalised Jordan form in the infinite dimensional case) and equivalent 
up to a more general (non necessary unitary) to the initial Hamiltonian. 
These Hamiltonians projected to the subspace S are Hermitian. 

If some of the eigenvalues are discrete and real, this same statement will 
obviously apply in the subspace of the corresponding eigenvectors. We shall 
study a simple example which illustrates this in a physical context, that of the 
Landau problem of a particle in two dimensions moving under the infuence 
of a constant magnetic field in the third direction [3] and show that this 
model exhibits the features which have been found for VT symmetric non- 
Hermitian systems. In particular, while the usual model is Hermitian, with 
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a real magnetic coupling, a deformation of the parameters in this model 
leads under certain systems to a non-Hermitian system with real eigenvalues 
and an imaginary coupling. In this respect it recalls to mind the work of 
Hollowood [1], who showed that affine Toda field theory with pure imaginary 
coupling possesses real energy levels. This system represents an anisotropic 
oscillator in interaction with a constant magnetic field. Some discussion is 
given of the case of transition between the Hermitian and non-Hermitian case, 
where two of the eigenvalues coincide, and the associated matrix cannot be 
fully diagonalised, but can be reduced by similarity transformations only to 
Jordan Normal form. 



2 Two-dimensional Case 

Take a 2 x 2 matrix 



M 



si s 3 
s 4 s 2 



(1) 



Then any finite dimensional matrix eigenvalue problem can always be trans- 
formed to an infinite dimensional Quantum Mechanical eigenvalue problem 
by the introduction of Fock space creation and annihilation operators aj, cij, 
so that matrix My turns into the Hamiltonian o,\Mija,j. This gives in 
this case the Hamiltonian 

H = s\a\ai + S2a|a2 + s^a\a2 + s±a\ai , (2) 

where the coefficients Sj, i = 1 . . . 4 are complex constants and two 
creation and annihilation operators satisfying the usual canonical commuta- 
tion relations, with only the following non zero commutators, 

didj — a^aj = 5ij . (3) 

The vacuum state | > is defined as usual as the normalized state which is 
annihilated by the a's 

a l \0>=O (4) 
where O is the zero norm state in the Hilbert space. 
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3 Diagonalisation 

Define 



A = ^/(s 1 -s 2 ) 2 + 4s 3 s 4 ; A± = ^(s 2 - si) ± ^A , (5) 
and construct the linear combinations 

at = ni [?~<A + 4^ (6) 

ct+ = n 2 \ P~a\ + a 2]"j (7) 

ai = ^(z ai+ A° 2 ) (8) 

1 /— s 4 — A_ \ , . 

a a = _^_ ax + - 5 -a a j . (9) 

In this construction, a+ denotes a creation operator, which should not be 
thought of as the Hermitian conjugate of atj which is also a destruction 
operator, since there is no Hermitian Hamiltonian (at ^ art"). However the 
properties 

< | Oia/ | > = 5ij (10) 

hold. Indeed the ay's and o^'s satisfy exactly the same commutation rela- 
tions as do the a's, as a consequence of the fact that the matrix (JTJ) can be 
diagonalised by a similarity transform. 

This implies that observables in the theory can be calculated, as n particle 
states will have the form 

I r, n - r > — 1 = (atY (atY " I > . (11) 
yjr\(n-r)\ V J K J 

These states are orthogonal to the following adjoint states 

| adj [r,n-r] >= (a{Y (4)""' I > (12) 

Jr\(n — r)! 

in the sense that 

< adj[p', g'] | p, g > = 5 p > p 5 q > q . (13) 
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The limitation of the use of non-Hermitian Hamiltonians is that the notion of 
+ conjugation and reality of eigenvalues is specific to the Hamiltonian used 
and is not universal as is the case with Hermitian conjugation, as Bender et 
al have remarked pQ. Hermiticity also guarantees reality of eigenvalues, inde- 
pendently of details of the Hamiltonian. This universality is a consequence of 
the fact that the inverse of a unitary matrix, which diagonalises a Hermitian 
matrix, is its own Hermitian conjugate, i.e. since IP = U^ 1 for a unitary 
matrix U the columns of are orthogonal to the rows of U, the operation 
of Hermitian conjugation works independently of the Hermitian matrix to be 
diagonalised. In the case of a 2 x 2 matrix H we have 



This construction shows that the n particle states have energies of the form 



where m runs from — n to n in steps of 2. There are several interesting 
features of this result. The eigenvalues will be real and distinct provided 
S1+S2 is real and (s±— S2) 2 +4s3S4 is positive, thus even when the Hamiltonian 
is non-Hermitian, the eigenvalues may be real. When this latter factor is 
zero, then the eigenvalues are degenerate. Then, either the Hamiltonian 
is proportional to the unit matrix (si = S2, S3 = S4 = 0) or, if at least 
one off-diagonal element is non zero, there is one eigenvector only and the 
Hamiltonian can be brought to the normal form (s[ = s' 2 , s' 3 = 0, s' A = 1). 
These cases are briefly discussed later. For non-Hermitian Hamiltonians, 
Bender et al have analysed the existence of a Parity operator V and a 
Conjugation matrix VT. These crucial issues will be taken up and generalised 
after we demonstrate a simple application to a physical example of the 2x2 
situation detailed above. 

4 Schrodinger interpretation 

The above problem is equivalent up to a constant energy shift to solving 
a Schrodinger equation for a particular Quantum Mechanical problem; set 



H 



Sia[ai + S2a 2 a2 + s^a^ + s^a^ai 





(15) 
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di = Pi + iqt etc; then the problem is equivalent to 

H - -»<£-»■> 

-(S3 + S 4)(^-^- - xy) + (s 3 - s A )(x— - y—) (16) 
axay ay ox 

up to a c-number addition. If S3 = — S4, then this is the Landau problem 

of a particle in a plane coupled to a magnetic field in the perpendicular 

direction, through its angular momentum with possibly an additional linear 

central force. To see this consider a constant magnetic field with potential 
B 

A = — ( x) in the Hamiltonian 



^(p'+p^ + (^) (z 2 + y 2 ) + ^-ip,x- Pl y) ] . (17) 



and make the identifications; s\ = S2 = l/(2m), and S3 = — S4 is pure imag- 
inary. This Hamiltonian is Hermitian with real eigenvalues. However, the 
choice of S3 = — S4 and real, eliminates the xy cross terms and also gives real 
eigenvalues, provided that the factor (si — s 2 ) 2 — 4s| is positive which nec- 
essarily entails that si 7^ s 2 . In this case the Hamiltonian is non-Hermitian 
and it represents a constant magnetic field coupled to an anisotropic oscil- 
lator with pure imaginary coupling. In all the examples Bender et al 2J[3 
have constructed with non-Hermitian Hamiltonians and real eigenvalues the 
Hamiltonian is VT symmetric. Dorey, Dunning and Tateo [Zj have recently 
given a proof that the spectra of a number of VT invariant Hamiltonians 
are entirely real. However, on its own VT invariance only implies real or 
complex conjugate pairs of eigenvalues. This is the situation under consider- 
ation here as VT invariance means invariance under (x, y) — > (— x, —y) and 
% — > —i. At the same time p does not change. The Hamiltonian will then be 
VT symmetric since B —> B, but the eigenvalues may be complex conjugate 
if the positivity condition is violated. 



5 Real Eigenvalue conditions 

As is well known, the eigenvalues of a Hermitian matrix are always real. In 
general the condition for the reality of eigenvalues depends on more specific 
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details of the matrix. However in the case of an n x n matrix H necessary 
conditions for the existence of real eigenvalues are given by 



Trace (if r ) = real; r = 1, 2, . . . , 



n . 



(18) 



This result follows simply from the observations that the characteristic poly- 
nomial of H, being an invariant under similarity transformations, has coef- 
ficients expressible in terms of the traces of powers of H and that if all the 
eigenvalues are real, must have the same eigenvalues as H . A further 
necessary and sufficient condition for an arbitrary 2x2 matrix to possess 
real eigenvalues is given by the requirement that 2Traceif 2 — (Traceif ) 2 > 0. 
This expression is just a translation into invariant form of the condition of 
section ©. Unfortunately, an invariant criterion even in the 3 dimensional 
case is rather complicated. When the traces are real and A = 0, the two 
dimensional matrix H may be expressed as 



where A is null or the zero matrix; i.e Traced = TraceA 2 = or A = 0. In 
the former case the matrix is equivalent by a change of basis to a Jordan 
Normal Form 



6 "Conjugation" and "Parity". Finite Di- 
mensional Case 

In this section, we extend the notions of "conjugation", and/or of "parity" 
introduced in [2] to arbitrary finite dimensional matrices representing the 
Hamiltonian with real eigenvalues. In certain cases the arguments can be 
extended immediately to infinite dimensional spaces (see Section (JHJ) for an 
example). It is also not difficult to extend them to the case where some of 
the eigenvalues are equal or are complex or when the starting matrix 
below has parts in a Jordan form 

We first recall a few well-known facts. Let Md be a real diagonal matrix in 
p dimensions with, for later simplicity, all elements different 




(19) 





(21) 
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and let N be an arbitrary invertible p x p complex matrix Let the scalar 
product of the column eigenvectors ip and be defined by the obvious 

< ijj | >= fltf, (22) 
where = (A*) tp is the Hermitian conjugate matrix. 

For any matrix M which has the eigenvalues \ii there exists a matrix N such 
that 

M .Y.U/.Y ; •< > M d — N~ x M N . (23) 
The p vectors <f>i, % — 1, . . . ,p, with components 

{<k)j = Sij , (24) 

are obviously eigenvectors of Md with eigenvalues /ij. They are orthonormal 
for the scalar product 

< (f)j | 4> k >= 8 jk . (25) 

The vectors tpi = Nfa are obviously eigenvectors of the general M with the 
same eigenvalue. 

We now present a few results providing, in the rather general case (j21)l . 
the general solution to the question of existence and construction of a suitable 
scalar product < <fi \ ip >a, of a "parity" operator P and of a "conjugation" 
operator C. 

1. The vectors ipi are orthonormal for the Hermitian scalar product (|26|) 
defined by 

<^|0> G =^t G ; with G = (N' 1 ^ N^ 1 . (26) 

Indeed 

Sjk = < 4>j I 4>k > 

= < iv-Vi I iv-Vfc > 

= < tpj | Vfc > G • (27) 

Because of the form of G (126(1 all the requirements of the scalar product 
are met. Note however that if iV is unitary, then G is the identity 



7 



matrix; otherwise, it depends upon the structure of N. This means that 
while Hermitian conjugation can be defined universally for Hermitian 
Hamiltonians. In the non-Hermitian case, it is a more specific matrix. 
Two matrices M lead to the same metric G if the iV which defines the 
first one is equal to the N which defines the second multiplied by an 
arbitrary unitary matrix. 

2. We denote complex by T the symbolic operator which performs the 
c-number conjugation, T 

TA = A*T. (28) 

There exist always a matrix P (generically called "parity" as Bender 
suggested) which has the following property 

(PT) M = M (PT) 

PM* = MP . (29) 

Take any invertible matrix K which commutes with M d , 

KM d = M d K . (30) 



In particular any diagonal matrix will do if the eigenvalues since are all 
different. If some eigenvalues are equal, K may belong to the stability 
group of M d . Then 

P = NK(N*y 1 . (31) 

Indeed 

PM* = MP 
P (N MdN^Y = N M d N~ x P by (j22D 

PN* M d (N~ 1 )* = N M d N~ l P by (j2U) 
(N~ 1 PN*)M d = M d (N- l PN*) . (32) 

Using (}3T]) and (jH2J), the statement is proved. 
3. The matrix C (conjugation) is defined as having the following property 

CM = MC 
C 2 = 1 . (33) 
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The general solution for C is 

C = NK.N- 1 (34) 

where K s is such it commutes with Md and K 2 = 1. In the generic 
case, K s is a diagonal matrix with elements ± only. Up to a trivial 
overall sign, there are 2 P_1 independent C's. Indeed the first equation 
of © gives 

CM = MC 
CNMaN' 1 = NM d N- l C by 
(N- l CN)M d = Md^N^CN) (35) 

which implies that N~ 1 C N = K commutes with M d - The second 
equation implies moreover that 

K 2 = I . (36) 

In the two dimensional example of section (j2J), C is either the identity 
matrix or the matrix 



st - s 2 2s 3 
2 S4 S2 — S\ 



(37) 



up to sign. 



7 More Bosonic Creation Operators 

In this section, we elaborate on two or more particle states obtained for the 
simplest Hamiltonian (J2J). We will treat in great detail the two and four 
particles states. 



7.1 Two Particle States 

Let us first elaborate on the two particle states constructed from the operators 
of section (©)• Suppose that we were to start with the products b\, i = 1, 2, 3 
of two creation operators 

b\ = (at) 2 
b\ = a\al 

bl = (4) 2 (38) 
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which construct the three two-particle states b\ | >, i = 1, 2, 3 when applied 
to the vacuum of the dj operators. The corresponding three-dimensional 
matrix M is 

/ 2 Sl s 3 \ 

2s 4 S! + s 2 2s 3 . (39) 
V s 4 2s 2 / 



It satisfies (with H from 



and hence 



if&l I 



M [2] & f 



Mjffo] | o > 



(40) 



(41) 



The diagonalisation of M provides the two particle spectrum. Using A from 
(JsjJ), we find, as expected, three eigenstates [3f which fulfill the eigenvalue 
equation (with H from (j2J)) 



H(3f | >=faPf I > 
They are up to a factor, obviously from (J2J) 



(42) 



-2s^t + 2s 3 (-A + si - s 2 )4 



(A(si - s 2 ) - (si - s 2 ) 2 - 2s 3 s 4 )4 



oc 



an 



-s 3 6| + (si + s 2 )& 2 + s 4 6. 



oc a^aj 



2sg&i + 2s 3 (A + si - s 2 )&2 



+ (-A( Sl - s 2 ) 
oc (a^) 2 



(si - s 2 ) 2 - 2s 3 s 4 )6 3 



(43) 



with eigenvalues (see (fT5|)) 



/'1 

/'3 



si + s 2 + A 

Sl + s 2 

si + s 2 - A 



E22 

E<2,0 
-E'2,-2 



2^i,i 

£a,i + #1,-1 
2^i,-i . 



(44) 



This discussion suggests the following question: how can this spectrum 
be obtained directly from the operators defining the two-particle space? We 
will first show that a naive approach has strong shortcomings. 
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1. A first naive answer is as follows. Introduce symbolically the new op- 
erators hi which are supposed, with the bp to satisfy the canonical 
commutation relations 

5 i3 (45) 



bi,b] 



[21 

and define a two particle Hamiltonian -ff naive as 

-'-'naive u i lvl ij u ] ■ 

Then obviously but, we insist, symbolically we obtain 

\H m 6 f ] = M [2l 6 f 

-'-'naive) i - LV± ii 



31 3 



and on a vacuum | 0^ > such that 

h I [2] >= 



we obtain 



ff£U I [2] >= M^b) i 0® > . 



(46) 



(47) 



(48) 



(49) 



It is not difficult to see that it is impossible to construct fej satisfying 
(}4l)j) out of products of two dj. However if the vacuum | [2] > is thought 
to be the a vacuum I > 



QPl >: = | > 



and the weaker condition 



b u fo]] | > = kb\ | > = 6u | > 



is imposed, the 6, can be identified to be 



b\ = -bx, bi = 62, bi — -63 in the weak sense 



(50) 



(51) 



(52) 



with, obviously, b\ = a\, b 2 = a\a 2 , b 3 = a\. The equation PH]) then 
holds. 

2. In a second approach, one tries to construct directly an quadratic 
both in a\ and in and which satisfies the basis equation (|4Tj) 



HM,b\] I >= M [ M I > 



(53) 
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as a result of the basic canonical commutation relations 

Such a Hamiltonian does not exist in the strong sense (when the vacuum 
is removed in (J53j0 . In the weak sense, it exists as 

H® = blNlfbj . (54) 

where the matrix is 





( 


«1 


S3 





N [2] = 




S 4 


Si + s 2 






\ 





s 4 


s 2 



This result is obviously equivalent to the naive approach once the 6j 
are expressed in terms of the bi (J52j) . 



7.2 Four Particle States 

Let us proceed in the same way for the four particle states which can be 
thought either directly as the four particle states in the original operators ( 

a\, i = l,2) m 

d\ = (a{) 5 - 1 (air 1 , i = l,...,5 (56) 

or as the compound states constructed out of two of the basic two particle 
states b], % = 1, 2, 3 



d\ = (b\f , d\ = b\bt , 4 = b\b\ 



d\ = (b\? , d\ = b\b\ , d\ = (4) 
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Let us note immediately that there are only five states built with the five 
d\ while there are six states built with the six d\. This is due to the obvious 
fact that 

d\ = d\ (58) 

when constructed from the a's. 

Let us treat in succession the d\ case and the d\ case. 

• The 5x5 matrix , analogous to (j3Uj) and defined by 

H,d{\ =Mf}d] , (59) 
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IS 



/ 4si 


S3 








\ 


4s 4 


3si + s 2 


2s 3 











3s 4 2 


>i + s 2 ) 


3s 3 











2s 4 


si + 3s 2 


4s 3 


V 








s 4 


4s 2 / 



(60) 



From this expression, it is easy to generalize to the general form for 
states with n original particles. Its eigenvalues and eigenstates, analo- 
gous to (jUIJ), (jUj) can be read off explicitly in (JTTJl and (fT5|). 

Note that, in the weak sense, there is a Hamiltonian ifW which may 
be written directly in terms of the cL and 



^ 4 UJ1 io>=mMio 



3* 



> 



and is 



with 



(61) 
(62) 



ATM 





Si 


S3 








\ 




s 4 


3si + s 2 


3s 3 













3s 4 


3(si + s 2 ) 


3s 3 













3 s 4 


si + 3s 2 


S 3 


I 











s 4 


s 2 y 



(63) 



Let us now try to define the 6x6 matrix in the same way for the 
dt 

121 Jtl 



naive; **« 



(64) 

again starting from the same canonical commutation relations for the 
&}'s and bj gSJ). We find 



M 



4sx 


s 3 














4s 4 


3si + s 2 


2s 3 


2s 3 











s 4 2 


>i + s 2 ) 





S3 








2s 4 





2(si + s 2 ) 


2s 3 











2s 4 


2s 4 


si + 3s 2 


4s 3 














s 4 


4s 2 



(65) 
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Since the d\ are perfectly defined there is no ambiguity. However, 
coherence with (|58|) implies that an arbitrary combination of 4 — d\ 
can be added to the right-hand side of (jMj) and Ad\f can be replaced 
by 

M\f = M®+r i 5 a -r j 8u t . (66) 

namely 



*? 














S3 














4s 4 


3 Si + s 2 


2 s 3 


2 s 3 








Ti 


«4 + r 2 


2 (si + s 2 ) + r 3 






r 6 


-n 


2 s 4 - r 2 


-r 3 


2(si + s 2 ) 


- r 4 2 s 3 - r 5 


-r 6 








2s 4 


2s 4 


Si + 3 s 2 


4 s 3 














s 4 


4s 2 



(67) 

Making the same combination for the left-hand side, it can easily be 
checked that 

41c 4 - 4] = (2si + 2s 2 + r 3 - r 4 ) (4 - 4) (68) 

i.e. a combination of 4 — 4 only. This result is coherent when 4 — 4 
is put to zero. 

The eigenvalues of the matrix M-j depend upon the values of the 
arbitrary parameters r; since the secular equation is 

(A-£ 4)4 ) (A-^4,-4) (A-^4,2) (A-^-2) (A-^4,0) (A-£ 4 , -r 3 +r 4 ) (69) 

When r 4 7^ r 3 , the eigenvalues are all different and the matrix can 
be diagonalised. When r 4 = r^, two eigenvalues become equal. Then 
there is always the eigenvector (0, 0, 1, —1, 0, 0) tp . The condition for the 
existence of a second eigenvector is 

r 6 s 4 + 2r 5 sis 4 - 2r 5 s 2 s 4 + r 3 (s 4 - s 2 ) 2 - 2r 3 s 3 s 4 

- 2r 2 sis 3 + 2r 2 s 2 s 3 + r^j = (70) 
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When this condition is satisfied, the matrix can be transformed to a 
fully diagonal form; otherwise it is reducible to a Jordan Normal form 



#4,4 



















#4,2 




















#4,0 


1 

















#4,0 




















#4,-2 




















#4, 



(71) 



8 Extension to a Fock space 

We shall now promote the discussion to the case of a Fock-space Hamiltonian 
linear in the products of one creation and one annihilation operator a\dj of 
bosons of two (or more) different species. The matrices N now become 
functions of the creation and annihilation operators (aj, al) which satisfy 
the usual canonical commutation relations (see ©)• 

In analogy with the 2-dimensional case ()14j) and with the discussion of 
the preceding section, the diagonal Hamiltonian is defined as 

H d = Yl Hj a ] a j ( 72 ) 
j 

and the related Hamiltonian H by 

H = NHaN- 1 (73) 

which should be linear in the products a\aj since this is a free field theory. 
Suppose 

NcnN- 1 = Y, ci S a i (=«i) 

3 

NajN- 1 = Y.^ a ) (=<4) • ( 74 ) 
j 

Indeed, we then find 
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(75) 



k rn 



If an interpretation in terms of creation and annihilation operators is to 
remain, the right hand sides of ()74|) should obey the same commutation rela- 
tions (jHJ) as the a's. This implies for the matrices c and d having respectively 
Cij and dij as components the restrictions 

d [i]t P = (c [i])-i (76) 

After some algebra, using (fTBj) and defining 

d® = - 1 , (77) 

one finds 



N,a 



df p N aj 
df} a) N . 



(78) 
(79) 



Using a power expansion of iV in terms of the a's, it is easy to show that the 
quadratic part of iV is 



Q = J2 d 



t 2 U 



and that the general N is constructed from Q as 

N =: exp(Q) : (81) 

where the symbol : . . . : denotes the normal product (annihilation operators 
written at the right of creation operators). Indeed, we have 

[di,N] = d : exp(Q) : 

a \ 
i 

= :exp(Q): (d f Q 



(80) 



Eq.(|7HJ) 



(82) 



but also 



N,aj 



d ai exp(Q) : 
(^Q)exp(Q) 



Y,d [ §a}]N <=; 



Eq.® 



(83) 
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This completes the proof. 

It is then not difficult to see that, inversely that any H of the form 



can be put into diagonal form Hd by using the inverse of N and adjusting 
the coefficient c^\d\ 1 ^ suitably. The /i's are the eigenvalues of the matrix 
constructed with the hij which is supposed to be diagonalisable at this stage. 

In particular for the case of p = 2, this can be read off directly for an 
arbitrary H from formulas © flSJ) © for the coefficient c [1] , of Eq.flS). 

Our arguments about the existence of a Hermitian Hamiltonian when the 
eigenvalues are real also applies to the infinite dimensional case. While we 
were concluding this article a paper by Mostafadazeh |Hj appeared, demon- 
strating that a quantum mechanical anharmonic oscillator with potential ix 3 
giving rise to a non-Hermitian Hamiltonian can be transformed to a Hermi- 
tian form, thus giving an alternative proof of the reality of energy eigenvalues. 



The Landau problem of a particle in a magnetic field has been shown to 
demonstrate the phenomenon of a non-Hermitian Hamiltonian giving rise to 
real eigenvalues when the coupling to the magnetic field becomes pure imagi- 
nary has been exhibited. All properties of this model are explicitly calculable 
as it is really a free field theory in disguise, which illuminates the conditions 
under which non-Hermitian Hamiltonians may have real eigenvalues. Indeed, 
starting from real eigenvalues in a finite system, which is then subject to a 
similarity transform, it has been shown how to define creation and annihi- 
lation operators, and their conjugation to build a Fock space for a general 
non-Hermitian Hamiltonian. In all the cases that we have studied in detail, 
we have observed that, at least formally, when a non-Hermitian Hamilto- 
nian possesses real eigenvalues, there exists a change of basis such that the 
transformed Hamiltonian in the corresponding subspace becomes Hermitian 
Other properties of such systems have been discussed, notably the conditions 
for degenerate eigenvalues. In the 2x2 case and extensions thereof this is 
a sign that the associated matrix is not fully diagonalisable and marks the 
transition between real and complex conjugate pairs of eigenvalues. 
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